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Abstract The first part of our analysis uses the wavelet me¬ 
thod to compare the Quantum Chromodynamic (QCD) pre¬ 
diction for the ratio of hadronic to muon cross sections in 
electron-positron collisions, R, with experimental data for R 
over a center of mass energy range up to about 7 GeV. A di¬ 
rect comparison of the raw experimental data and the QCD 
prediction is difficult because the data have a wide range of 
structures and large statistical errors and the QCD descrip¬ 
tion contains sharp quark-antiquark thresholds. However, a 
meaningful comparison can be made if a type of “smearing” 
procedure is used to smooth out rapid variations in both the 
theoretical and experimental values of R. A wavelet analysis 
(WA) can be used to achieve this smearing effect. The sec¬ 
ond part of the analysis concentrates on the 3.0 - 6.0 GeV 
energy region which includes the relatively wide charmo¬ 
nium resonances y/( 1“). We use the wavelet methodology to 
distinguish these resonances from experimental noise, back¬ 
ground and from each other, allowing a reliable determina¬ 
tion of the parameters of these states. Both analyses are ex¬ 
amples of the usefulness of WA in extracting information in 
a model independent way from high energy physics data. 


1 Wavelet Transformations 

Let us start with a brief description of the continuous wave¬ 
let transformation (WT). The WT of function /(f) is defined 
by 

w(a,t) = -2= / (p* (- - ^\f(t')dt', (1) 

yjaCq, J \ a J 

where C l(l is a normalization constant subject to the choice 
of wavelet. The decomposition described by equation (1) is 
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performed by convolution of the function /(f) with a bi- 
parametric family of self-similar functions generated by di¬ 
latation and translation of the analyzing function <p(f) called 
a wavelet, 

= V , ( 2 ) 

where the scale parameter a characterizes the dilatation, and 
b characterizes the translation. It is a kind of “window func¬ 
tion” with a non-constant window width. High frequency 
wavelets are narrow due to the factor 1/a, while low fre¬ 
quency wavelets are much broader. The function <p(f) should 
be well localized in both time and Fourier space and must 

+°o 

obey the admissibility condition, f (p(t)dt. This condition 

requires <p(f) must be an oscillatory function and, if the in¬ 
tegral (1) converges, the completeness of the wavelet func¬ 
tions provides the existence of inverse transformation. 
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In contrast to Fourier analysis, the WT depends both on 
f and the frequency providing an optimal compromise with 
the uncertainty principal. One of the advantages of wavelet 
analysis is a fairly low sensitivity of the restored signal to 
any physically reasonable continuation of the function /(f) 
outside the interval (t m j n ,tmax) where the data are known. 
To fill in gaps between the experimental points we use a 
linear interpolation (different interpolations lead to minimal 
difference in the restored signal). Note that since the average 
value of any wavelet is zero, the mean value of the WT is 
zero, so that (/) must be added to the reconstructed signal 
to restore the mean value of the original signal. 

Wavelets with good localization and a small number of 
oscillations are commonly used to recognize the local fea¬ 
tures of data, and to find the parameters of dominating struc¬ 
tures (location and scale/width). In this work, we use one of 
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the most popular wavelets of this type, the so-called “Mexi¬ 
can Hat”, 

(p(t)=(l-t 2 )e-‘^ 2 . (4) 

This wavelet is plotted in Figure 1 for three values of the 
scale parameter, a = 1, 2 and 0.5. 



Fig. 1 Wavelet “Mexican Hat” q>(t) = (1— t 2 ) e 'I 2 . Solid line — 
<p(f), dashed line — (p(t/ 2), dotted line — <p(t/ 0.5). 


2 Wavelets and the R ratio 

Our first goal is to use wavelet methodology to compare the 
QCD prediction of R to experimental data in the center of 
mass energy range up to about 7 GeV. In zeroth order in the 
strong coupling constant a s , the ratio R is given by 


m) 


o (e + e hadrons) 

a(e + e~ —> ) 


3^e 2 q =R (0) (Q), 


(5) 


where the summation extends over the quark flavors q = 
u,d,s,c available up to the center of mass energy Q. The 
QCD a s corrections to r(°\Q) can be presented in different 
forms. To be specific we adopt the form of reference [1], 

R(Q)=3^e 2 q T(v q )[l+g(v q )^], ( 6 ) 

<3 


where 


_ a s 
3 ?= — 

n 


1 + —+C 2 (— +C 3 (— +... 

n \ n / V n ) 

v? = [! - 4 m 2 q /Q 2 } 1/2 , T(v) = v (3 - v 2 ) /2 


4 K 

’ n 

3 + v 

( 71 

Ml 

3 

_ 2 v 

4 

V 2 

4 K J 


and the summation in ( 6 ) is over all quark flavors whose 
masses m q are less than Q/2. For this analysis we use the 
quark masses, the form of the coefficients C 2 and C 3 calcu¬ 
lated in [3] and the energy dependence of a s (Q 2 ) from [ 2 ]. 

It is very hard to determine the role of these QCD cor¬ 
rections by comparison with experimental data for R(Q) due 
to the large statistical errors and plethora of overlapping and 
interfering resonances in this region. An additional compli¬ 
cating factor is that the QCD perturbative approach exhibits 


sharp quark-antiquark thresholds. However, a meaningful 
comparison can be made by applying some type of “smear¬ 
ing” procedure, which has the effect of smoothing out rapid 
variations in both the theoretical and experimental values of 
R. 

Before describing our wavelet analysis of this data it is 
worth considering the methodology developed in references 
[1] and [2] to compare experimental data and QCD predic¬ 
tions. The smearing procedure used in these analyses calcu¬ 
lates a smeared ratio R as follows, 


R(s, A) 


A 

K 


Smax , 

f 

7 (V — s ) 2 + A 2 


(7) 


where \fs = Q is the square of the center of mass energy and 
A is the “smearing” parameter. In references [1] and [2] to 
evaluate the integral (7), it was necessary to exclude sharp 
resonances, such as the i/(3.100) and the much wider p 
peak. In addition, a term is added to account for the contribu¬ 
tion from s max to °°, assuming that R remains constant above 
Smax ~ 60 GeV 2 . Originally, the smearing procedure in ref¬ 
erence [1] supposes a global constant value A = 3 GeV 2 in 
(7). However it was found in reference [2] that for different 
energy regions it would be better to use different values of A. 
Note that the use of an energy dependent A in (7) reflects the 
necessity for different treatment of different energy scales. 

The WT methodology provides an alternative, model in¬ 
dependent, smearing method, which does not require differ¬ 
ent treatment in differing energy scales. Under WT to sep¬ 
arate the signal from the background noise, wavelet recon¬ 
struction is performed for scales greater than a certain scale 
a noise ~ the boundary, or cut-off, scale [4]. In deciding on the 
appropriate boundary scale that will separate the noise-like 
high frequency components of the data we take a pragmatic 
line of reasoning. That is, the best choice for a no ise is the 
smallest value which will smooth out any rapid variations in 
the data enabling us to reproduce stable results for low fre¬ 
quencies (resonance area). A similar pragmatic strategy was 
applied in the analyses of [ 1 ] and [ 2 ] in choosing the param¬ 
eter A of equation (7). The best value of A is large enough to 
compare the smeared R with QCD models, but not so large 
that all the fine detail of the data is smoothed away. 

Figure 2 displays wavelet reconstructed (smeared) ex¬ 
perimental data with cut-off values a no i se = 0.6 (bold dashed 
curve) and a no i se = 1.2 (bold solid curve). The data is a com¬ 
pilation of measurements from many different experiments 
obtained from the Particle Data Group (PDG) [5]. In con¬ 
trast to the analyses of [1] and [2], under the WT approach 
there is no need to remove sharp resonances (1/(3.100) etc.) 
by hand, they effectively become part of the high frequency 
noise and the wavelet analysis smears them out. The thin 
dashed line and thin solid line are wavelet reconstructed 
QCD curves with a no i se = 0.6 and a no i se = 1.2, respectively. 
As can be seen, the two curves with a no i se = 1.2, repre- 
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Fig. 2 Wavelet reconstruction of R e + e -. Dashed and solid bold curves 
correspond to wavelet reconstruction of the experimental data with 
a noise = 0.6 and a no i se = 1.2, respectively. Dashed and solid thin curves 
correspond to wavelet reconstruction of QCD calculations (6) with the 
same a„ 0 , je values. 

senting experiment and theory, are in good agreement. It 
should be noted that the contribution to all of the restored 
data curves in Figure 2 from the data region above 6.5 GeV 
is negligible. 


3 y/ states above the DD threshold and the wavelet 
procedure 

Detailed information on the charmonium resonances t/r(l _ ) 
above the DD threshold at 3.73 GeV comes primarily from 
the measurement of R e + e -. This data, provided by the Par¬ 
ticle Physics Data Group (PDG) [5], comes from the work 
of many experimental collaborations over a period of more 
than 30 years. Several broad vector resonances are observed 
with varying degrees of clarity. The current best estimate of 
the masses and widths from the PDG are presented in Table 
1 . 

The main difficulty encountered by each of the colla¬ 
borations in making these measurements is large statisti¬ 
cal errors and hence the difficulty of separating resonance 
“signal” from noise. There is significant disagreement be¬ 
tween the collaborations on many of the resonance parame¬ 
ters. Therefore, due to the possibility of systematic errors we 
do not combine measurements from different experiments, 
but choose to base this initial analysis on data from the BES 
collaboration which exhibits clear evidence of four broad 
resonances above the DD threshold. In addition to the PDG 
values, Table 1 shows the masses and widths of these reso¬ 
nances from the most recent BES analysis [6]. 

Before analyzing the charmonium resonances parame¬ 
ters, we first determine the non-resonant background con¬ 
tribution using a WT followed by wavelet reconstruction. 
This methodology, with its excellent scaling property, al¬ 
lowing the analysis of data with varying resolution, is ide¬ 
ally suited to separate resonances from noise, background 
and each other. A very helpful representation of the WT 
revealing all the features of the complete spectrum of the 


signal is the “time- frequency” plane (Figure 3(a)). This is 
a multi-resolution spectrogram, which shows the frequency 
(scale) contents of the signal as a function of energy. Each 
pixel on the spectrogram represents w(a.t) for a particular 
a (scale) and t (in our case t is energy, E = Q). The loca¬ 
tion of spots on the vertical axis (scale axis, a) corresponds 
to the width of the maximum. The intensity of dark spots 
shows the amplitudes of maxima. The WT image (wavelet 
plane) of the BES charmonium data obtained with the “Mex¬ 
ican Hat” wavelet is shown in panel (a) of Figure 3. The 
WT localizes the structures in a fashion that allows us to 
estimate the masses of the resonances and their widths - all 
four y/ resonances are clearly seen on the wavelet plane. The 
straight horizontal line corresponds to the boundary scale 
a no ise , which can be chosen to cut off small scale structures, 
which in this case corresponds to experimental noise. By 
choosing a larger value of a no ise it is possible to also cut 
out the resonances, leaving only the background. The back¬ 
ground curves for three substantially different choices of 
anoise (0.15, 0.25 and 0.35) are displayed in Figure 3(b). In 
the following analysis we show that the charmonium reso¬ 
nance parameters are not sensitive to the choice of a no j se . 



(*) 



Fig. 3 (a) Wavelet plane for a 6 [0.01. 1]; (b) background for a no j se = 
0.15 (solid line), a no i se = 0.25 (dashed line), a no i se = 0.35 (dotted line). 

For each value of a no i se we subtract this “wavelet back¬ 
ground” from the raw experimental data leaving only the 
contribution from the broad charmonium resonances. We 
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Table 1 Properties of vector (charmonium) resonances above the charm threshold from the PDG and the BES [6] collaboration 


Resonance 

Mass (GeV) 

Full Width (GeV) 

PDG [5] | BES [6] 

PDG 151 | BES [6] 

yr(3.770) 

3.773± 0.0003 3.772 ± 0.002 

0.027 ±0.001 | 0.030 ±0.009 

yr(4.040) 

4.039 ±0.001 | 4.040 ±0.004 

0.080 ±0.010 | 0.085 ±0.012 

yr(4.160) 

4.191 ±0.005 4.192 ±0.007 

0.103 ± 0.008 0.072 ±0.012 

yr(4.415) 

4.421 ±0.004 | 4.415 ±0.008 

0.062 ±0.020 | 0.072 ±0.019 


then perform a least squares fit of the “signal” in this en¬ 
ergy range to the sum of four Breit-Wigner resonances of 
the form. 


Rn 


— 2 


2 t _, 2 


Mir,. 


'em r — 1 


{s-M?y+M?r r 2 


( 8 ) 


where M r , 77, B/, and B/,,. are, respectively, the mass, total 
width, leptonic branching fraction and hadronic branching 
fraction of the resonance. Each Breit-Wigner has three fitted 
parameters, M, F and the product of 77/,. and B/ lr . It is worth 
noting that since these resonances partially overlap, if their 
decay channels are specified, we can improve their resolu¬ 
tion by using the multi-channel unitary scheme described in 
reference [7]. We defer this investigation to a future analysis. 

The values of the fitted parameters for each of the four 
Breit-Wigner resonances for backgrounds obtained with three 
different values of ci, WLU , are presented in Tables 2-4. The chi 
squared per degree of freedom for the three fits are 1.11, 


Table 4 Breit-Wigner fitted parameters for a no i se = 0.35 


Resonance 

Mass (GeV) 

Width (GeV) 

B,B h - 10 5 

V/(3.770) 

3.773 ± 0.0001 

0.022 ± 0.004 

0.888 ± 0.090 

yr(4.040) 

4.043 ± 0.003 

0.097 ±0.011 

1.155 ±0.061 

V/(4.160) 

4.165 ±0.004 

0.083 ± 0.013 

0.907 ± 0.073 

yr(4.415) 

4.423 ± 0.004 

0.090 ±0.014 

0.889 ±0.071 


ent for the three values of a no i se . Therefore, in Figure 4 we 
present only the fitted curve (8) for a no i se = 0.35. 

It should be noted that in this analysis, after subtracting 
the wavelet background, we exclude the four highest energy 
data points from the fitting procedure. Without these points 
the i//(4.415) resonance is “well-shaped”; if these four points 
are included the resonance is no longer “well-shaped” and 
its width is about twice that of the values presented in Table 
1 . 


1.22 and 0.99 for the a no is e values of 0.15, 0.25 and 0.35, 
respectively. 


Table 2 Breit-Wigner fitted parameters for a, w i se = 0.15 


Resonance 

Mass (GeV) 

Width (GeV) 

B,B h ■ 10 5 

V/(3.770) 

3.772 ± 0.0002 

0.032 ± 0.005 

0.983 ± 0.076 

V/(4.040) 

4.042 ± 0.001 

0.088 ± 0.013 

0.921 ± 0.067 

y/( 4.160) 

4.161 ±0.006 

0.100 ± 0.019 

0.664 ± 0.069 

V/(4.415) 

4.430 ± 0.007 

0.098 ± 0.014 

0.936 ± 0.071 



Fig. 4 Breit-Wigner fit of BES data with “wavelet background” sub¬ 
tracted, a no i se = 0.35 


Table 3 Breit-Wigner fitted parameters for a, w i se = 0.25 


Resonance 

Mass (GeV) 

Width (GeV) 

B,B h ■ 10 5 

V/(3.770) 

3.772 ± 0.0001 

0.022 ± 0.004 

0.888 ± 0.094 

V/(4.040) 

4.043 ± 0.003 

0.097 ±0.011 

1.110 ± 0.064 

V/(4.160) 

4.185 ±0.005 

0.083 ±0.016 

0.715 ± 0.077 

V/(4.415) 

4.423 ± 0.004 

0.090 ± 0.016 

0.819 ± 0.076 


As can be seen in Tables 2-4, the fitted parameters of 
the four charmonium resonances are not significantly differ- 


Justification for excluding the four highest energy points 
can be made as follows. If instead of the wavelet fitted back¬ 
ground we fit the data with four Breit-Wigners and a parabolic 
(or linear) background, including the four highest energy 
data points, we see that the fourth resonance is not “well¬ 
shaped” and the last few data points appear to be associated 
with the background rather than the 4th resonance, see Fig¬ 
ure 5. Furthermore, the fitted width of this fourth resonance 
(Table 5) is significantly larger than that of the second and 
third resonances and the published PDG value. In order to 
show the association of the highest energy points with the 
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parabolic background we do not subtract the background in 
Figure 5. It should be emphasized that inclusion or exclu¬ 
sion of these four data points does not significantly alter the 
background curve obtained via the wavelet method and that 
the choice of a polynomial background is arbitrary, whereas 
the wavelet background is obtained from the data itself in a 
model independent manner. 


Table 5 Breit-Wigner parameters for fit with parabolic background 


Resonance 

Mass (GeV) 

Width (GeV) 

B,B h ■ 10 s 

1/(3.770) 

1/(4.040) 

1/(4.160) 

1/(4.415) 

3.770 ± 0.001 
4.048 ± 0.003 
4.166 ± 0.003 
4.426 ± 0.005 

0.021 ± 0.004 
0.107 + 0.010 
0.086 ± 0.013 
0.148 + 0.022 

1.202 + 0.097 
1.946 + 0.079 
1.335 ± 0.067 
1.001 ± 0.085 



Fig. 5 Breit-Wigner and parabolic background fit of BES data (con¬ 
trary to Fig.4 the background is not subtracted from the data). 

To demonstrate the reliability of our WT methodology in 
extracting the parameters of these charmonium resonances 
we performed a study using simulated data as follows. Start¬ 
ing from the PDG values of masses and widths of the four 
charmonium resonances given in Table 1 we generate a spec¬ 
trum of four Breit-Wigners of the basic form of equation 
(8). To simulate real experimental data we discretize the 
spectrum from 2.0 to 5.5 GeV in 0.05 GeV intervals, then 
add a random vertical shift from -0.2 to +0.2 and random 
error bars in the range 0.15 to 0.4 to the data points. To 
complete the simulated data spectrum a representative back¬ 
ground curve is added to the data points. Finally we pass this 
simulated data through the same WT analysis chain applied 
to the real data. This gives us a wavelet background curve 
and the masses and widths of the Breit-Wigner resonances 
after the wavelet background subtracted fit. In addition to 
applying our WT methodology to the simulated data we per¬ 
form a Breit-Wigner and parabolic background fit in exactly 
the same way as we did for the real data. 

This study was performed for several different initial 
background curves. In all such cases the WT methodology 
accurately reproduced the input background shape. Fitting 


the wavelet background subtracted simulated data, as de¬ 
scribed above, also accurately reproduces the masses and 
widths of the assigned Breit-Wigner resonances, indepen¬ 
dent of the particular input background curve. The parabolic 
background least squares fit of the simulated data was also 
able to reproduce the Breit-Wigner masses and widths, but 
in contrast led to fitted backgrounds very different from the 
input background. We believe the fact that the masses and 
widths were accurately reproduced, despite the fitted back¬ 
ground being very different from the input background, is 
largely due to the clean nature of the simulated data. Real 
experimental data is clearly more complex, in which case the 
inability of the standard parabolic background least squares 
fit to extract an accurate background could lead to less reli¬ 
able resonance parameters. 


4 Conclusion 

Experimental measurements of the ratio, R, comprise a wide 
range of structures with large statistical errors, making di¬ 
rect comparison with the predictions of QCD very difficult. 
However, a meaningful comparison can be made provided 
that some kind of “smearing” procedure, similar to that de¬ 
scribed in [1], is used to smooth out rapid variations in R. 
A wavelet analysis can be used to achieve this smearing 
effect. We compare the WT of the predictions of perturba¬ 
tive QCD and experimental R data. The wavelet reconstruc¬ 
tion of the R experimental data preserves its main features, 
but with damped statistical errors and threshold singulari¬ 
ties. The WT of QCD perturbation theory is in good general 
agreement with the WT experimental data. 

Using the wavelet methodology to obtain the backgro¬ 
und in the charmonium energy range above the naked charm 
threshold provides an important, model independent, alter¬ 
native to other accepted methods. The masses and widths of 
the four vector mesons above the charm threshold, 1/(3.700), 
t//(4.040), 1/(4.160) and 1/(4.415), obtained from fitting wa¬ 
velet background subtracted data from the BES experiment, 
are found to be largely insensitive to the specific choice of 
the WT parameters and consistent with the BES and PDG 
reported values. 
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